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Abstract Many quasars and active galactic nuclei 
(AGN) appear in radio, optical, and X-ray maps, as 
a bright nuclear sources from which emerge single or 
double long, thin jets. When observed with high an- 
gular resolution these jets show structure with bright 
knots separated by relatively dark regions. Nonther- 
mal nature of a jet radiation is well explained as the 
synchrotron radiation of the relativistic electrons in an 
ordered magnetic field. We consider magnetic collima- 
tion, connected with torsional oscillations of a cylin- 
der with elongated magnetic field, and periodically dis- 
tributed initial rotation around the cylinder axis. The 
stabilizing azimuthal magnetic field is created here by 
torsional oscillations, where charge separation is not 
necessary. Approximate simplified model is developed. 
Ordinary differential equation is derived, and solved nu- 
merically, what gives a possibility to estimate quanti- 
tatively the range of parameters where jets may be sta- 
bilized by torsional oscillations. 



1 Introduction 

Objects of different scale and nature in the universe: 
from young and very old stars to galactic nuclei show 
existence of collimated outbursts - jets. Geometrical 
sizes of jets lay between parsecs and megaparsecs. The 
origin of jets is not well understood and only several 
qualitative mechanisms are proposed which are not jus- 
tified by calculations. Theory of jets must give answers 
to three main questions: how jets are formed? how 
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are they stabilized? how do they radiate? The last 
question is related to the problem of the origin if rel- 
ativistic particles in outbursts from AGN, where syn- 
chrotron emission is observed. Relativistic particles, 
ejected from the central machine rapidly loose their en- 
ergy so the problem arises of particle acceleration inside 
the jet, see review of Bisnovatyi-Kogan (1993). 

It is convenient sometimes to investigate jets in a 
simple model of infinitely long circular cylinder, Chan- 
drasekhar & Fermi (1953). The magnetic field in the 
collimated jets determines its direction, and the axial 
current may stabilize the jet's elongated form at large 
distances from the source (e.g. in AGNs), (Bisnovatyi- 
Kogan et al., 1969). When observed with high angu- 
lar resolution these jets show a structure with bright 
knots separated by relatively dark regions (Thomson 
et al., 1993). High percentages of polarization, some- 
times exceeding 50%, indicate the nonthermal nature 
of the radiation, explained as synchrotron emission of 
the relativistic electrons in an ordered magnetic field. 

Here we consider stabilization of a jet by magne- 
tohydrodynamic mechanism, connected with torsional 
oscillations. We suggest that the matter in the jet is 
rotating, and different parts of the jet rotate in dif- 
ferent directions. Such distribution of the rotational 
velocity produces azimuthal magnetic field, which pre- 
vents a disruption of the jet. The jet remains to be in a 
dynamical equilibrium, when it is representing a peri- 
odical, or quasiperiodical structure along the axis, and 
its radius is oscillating with time all along the axis. The 
space and time period of oscillations depend on the con- 
ditions at jet formation: the length scale, the amplitude 
of the rotational velocity, and the strength of the mag- 
netic field. The time period of oscillations should be 
obtained during construction of the dynamical model, 
what also should show at which input parameters may 
exist a long jet, stabilized by torsional oscillations. 
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Here we construct a simplified model of this phenom- 
ena, which confirms the reality of such stabilization, see 
also Bisnovatyi-Kogan (2007). 



the cylinder axis (z), and we approximate density by a 
function p(t, z), suggesting a uniform density along the 
radius. The components of the velocity and magnetic 
field are approximated as 



2 General picture 

Consider a long cylinder with a magnetic field directed 
along its axis (Fig.l). This cylinder will expand unlim- 
itly under the action of pressure and magnetic forces. 
The limitation of the radius of this cylinder could be 
possible in dynamic state, when the whole cylinder 
undergoes magneto-torsional oscillations. Such oscil- 
lations produce toroidal field, which prevent a radial 
expansion. There is a competition between the in- 
duced toroidal field, compressing the cylinder in ra- 
dial direction, and gas pressure, together with the field 
along the cylinder axis (poloidal), tending to increase 
its radius. During magneto-torsional oscillations there 
are phases, when either compression or expansion forces 
prevail, and, depending on the input parameters, we 
may expect tree kinds of a behavior of such cylinder. 

1. The oscillation amplitude is low, so the cylinder 
suffers unlimited expansion (no confinement) 

2. The oscillation amplitude is too high, so the pinch 
action of the toroidal field destroys the cylinder, and 
leads to formation of separated blobs. 

3. The oscillation amplitude is moderate, so the 
cylinder survives for an unlimited time, and its parame- 
ters (radius, density, magnetic field etc.) change period- 
ically, or quasi-periodically in time. Here we try to find 
a simple approximate way for obtaining a qualitative 
answer, and to make a rough estimation of parameters 
leading to different regimes. 

3 Profiling in axially symmetric MHD 
equations 



v r = ra(t,z), v v =rQ,{t,z), v z = 0; (I) 

B r = rh r (t,z), B ip = rh ip {t,z),B z = B z {t,z). (2) 
In this case the current components are written as 

cr dh v . cr dh r . ch v 

After neglecting velocity v z along the axis, we should 
omit the corresponding Euler equation, and the radial 
pressure gradient is approximated by the linear func- 
tion, what gives, when using appropriate dimensions 

where the constant A ~ I is connected with the equa- 
tion of state, P(t, z) is the pressure, R(t, z) is the ra- 
dius of the cylinder. In the subsequent consideration we 
consider an adiabatic case with a polytropic equation 
of state P — K p 1 . Neglecting the z derivatives in the 
Poisson equation, we obtain ipo = nGpr 2 . Substituting 
CH)— <j4j) into the original system of the equations, we ob- 
tain for the profiling functions the following equations 

da 9 , P „ _ 1 dh r „, 9 \ 

m +a - n = v - 2nGp+ *tp r ^ - 2hl ) ' 

(5) 



We simplify the system of MHD equations with axial 
symmetry, g^=0, for the perfect gas at infinite conduc- 
tivity, written in cylindric coordinates (r, <f>, z) (Landau 
and Lifshits, 1982). We use for this purpose a profil- 
ing procedure. There is no gravity in the direction of 
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(6) 
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d(aB z ) 
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do dR 

Fig. 1. — Jet confinement by magneto-torsional oscilla- — = —2ap, = aR. 

tions (qualitative picture) (>i 



(8) 
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It follows from ([H]),© relations, representing conserva- 
tion of mass, and magnetic flux equivalent to freezing 
condition 



dh v _ fd{Q/R 2 ) 

~dT = Cb 



dz 



z=0 



dR 

- 2ah vt — = aR, (14) 



P R 2 = C m (z), B Z R 2 = C b (z), B z = ^\p. 



(10) 



pR 2 = C m , B z R 2 = C b , B z = Q-p. (15) 
Initial conditions for the system (jTTjl - (TT5")) are 



In our subsequent consideration the arbitrary functions 
will be taken as constants: C m (z) = C m , Cb(z) = Ch- 
it was shown by Bisnovatyi-Kogan (2007)that the 
approximate system of equations describes correctly 
small perturbations, connected with radial and tor- 
sional modes, in the equilibrium self-gravitating cylin- 
der. We expect therefore that these modes will be de- 
scribed correctly in general nonlinear case in the jets 
where gravity is neglected. 



R — Ro, p — po 



Cm 
R l 1 



B, 



(16) 



5 = h r = h v = at t = 0. 
In (fTT|) - (fTS"]) we have used relations 



R a - C b R 



(17) 



4 Farther simplification: reducing the problem 
to ordinary differential equation 

For the relativistic jet we neglect gravity, without which 
the equilibrium static state of the cylinder does not ex- 
ist. We need to solve numerically the system of nonlin- 
ear equations (f5|)- (|10[) to check the possibility of the ex- 
istence of a cylinder, which radius remains finite due to 
torsional oscillations. Instead we reduce the system to 
ordinary equations. Supposing a constant bulk motion 
velocity along z-axis, we consider axially symmetric jet 
in the comoving coordinate frame. If the confinement is 
reached due to standing magneto-torsional oscillations, 
there are points along z-axis where rotational velocity 
remains zero in this frame. Taking Q = in the plane 
z = 0, let us consider standing wave torsional oscilla- 
tions with the space period zq along z axis. Then nodes 
with O = are situated at z = ±n^, n = 0, 1,2, .... 
Let us write the equations, describing the cylinder be- 
havior in the plane z = 0, where Q = 0. All values in 
this plane we denote by ("). We take also for simplicity 
A = 1. We have than equations in the plane z = as 



da n K 
a 2 = — 



C h fdh 



dt 



R2 AnC m V dz 



z=0 



h%R 2 



(11) 



valid for any time. If the cylinder rotational velocity is 
antisymmetric relative to the plane z = 0, Q = 0, and 
cylinder density distribution is symmetric relative to 
this plane, then we have extremum (maximum) of the 

azimuthal magnetic field h$, with (J^f^j ^ — 0, and 

zero value of h r = 0, which reaches an extremum (min- 
imum) in this plane with (^p) z _ = 0. The product 
ap also reaches an extremum in the plane z = 0, so that 
= 0. The term with z derivative in the 



z=0 



equation Q14p is not equal to zero, and changes period- 
ically during the torsional oscillations. We substitute 
approximately the derivative d/dz by the ratio 1/zo, 
where zq is the space period of the torsional oscillations 
along z axis. While £1 = in the plane z = 0, its deriva- 
tive along z is changing periodically with an amplitude 
f^o, and frequency uj, which should be found from the 
solution of the problem. We approximate therefore 



d{n/R 2 



= — COSldt. 

z R 2 



(18) 



Finally, we have from (fTT |) . (ll4p .(fl~5 |) the following ap- 
proximate system of equations, describing the non- 
linear torsional oscillations of the cylinder at given Zq 
and £Iq. 



C b fdh 



dz 



dK_ c ( d{a/R 2 ) 
dt \ dz 



h v h r R 2 
2nC m 



z=0 



0. 



(12) 
(13) 



da 
~dt 



a 2 



K 
R 2 



K R2 

27rC m 



dh v Q ~ 

— r^- = Ch — — cos ujz — Zan 
dt Zo R2 



dR 

~dt 



= dR. 



(19) 
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Fig. 4. — Time dependence of non-dimensional radius 
y (upper curve), and non-dimensional velocity z (lower 
curve), for D = 2.15 during a long time period. 




y (upper curve), and non-dimensional velocity z (lower 
curve), for D = 2.55. 



Fig. 10. — Time dependence of non-dimensional radius 
y (upper curve), and non-dimensional velocity z (lower 
curve), for D = 3.1, during a long time period. 
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6 Discussion 



Consider isothermal equation of state P = Kp = v 2 p, 
Vg < c 2 /3, v s is the sound speed. For ultrarelativistic 
pair-plasma K — c 2 /3. The parameter D, as a function 
of Rq, zq, initial po and B z q, J7q and is written as 



D = 



1 



r2 p2(-)2 



27rpo 2gW 2 t> 2 



(23) 



The amplitude of oscillations f2, and lo should be found 
from the solution of the nonlinear system ([5|)- (TlO|) . to- 
gether with the interval of values of D at which confine- 
ment happens. In the approximate system (f2"2"|) only 
D characterizes different regimes, what for given i?o, 
zo, po, and B z q determines a function Qo(D,lo) for the 
collimated jet. To find approximately a self-consistent 
model with D,q,lo(D) we use the the frequency linear 
oscillations lo = IzVa- The frequency of non- linear os- 
cillations is smaller, so we may write 
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= at k 2 Vi = 



Po4 



a n < 1. 



(24) 



therefore fi 2 i?§ = 2tt 2 Da 2 n v 2 s < c 2 , R$ = 



K 



.,2 PoK 



o a2 ^ g 2 . On the edge of the cylinder the rotational 
velocity cannot exceed the light velocity, so the solu- 
tion has a physical sense only at v 2 < 27r 2 Da 2 ~ 4^2"- 
Taking a 2 = 0.1 for non- linear oscillations we obtain 
a restriction v^ < To have the sound velocity not 
exceeding c/2, the jet should contain baryons, which 
density po exceeding about 30% of the total density of 
the jet. In a dense quasispherical stellar cluster around 
a supermassive black hole, the accretion disk is chang- 
ing its direction of rotation, due to different sign of the 
angular momentum of the falling stars. In this sutu- 
ation the magneto-torsional oscillations should be in- 
evitably generated in the outflowing jets. The knots 
in jets are observed in different objects, including a fa- 
mous M87 jet. In the last one the knots are visible in 
radio, optics (Hiltner, 1959), and X-rays (Wilson and 
Yang, 2002). It is possible that they are connected with 
magneto-torsional oscillations. 
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